Quantitative analysis of pedestrian counterflow in a cellular automaton model 
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Pedestrian dynamics exhibits various collective phenomena. Here we study bidirectional pedes- 
trian flow in a floor field cellular automaton model. Under certain conditions, lane formation is 
observed. Although it has often been studied qualitatively, e.g. as a test for the realism of a model, 
there are almost no quantitative results, neither empirically nor theoretically. As basis for a quan- 
titative analysis we introduce an order parameter which is adopted from the analysis of colloidal 
suspensions. This allows to determine a phase diagram for the system where four different states 
(free flow, disorder, lanes, gridlock) can be distinguished. Although the number of lanes formed is 
fluctuating, lanes are characterized by a typical density. It is found that the basic floor field model 
overestimates the tendency towards a gridlock compared to experimental bounds. Therefore an 
anticipation mechanism is introduced which reduces the jamming probability. 



I. INTRODUCTION 

From a physical point of view, human crowds are 
many-body systems which show many interesting col- 
lective effects. Examples are jamming at bottlenecks, 
density waves, flow oscillations, pattern- and lane forma- 
tion [l| . The latter terms the effect that pedestrians in a 
crowd segregate spontaneously according to their desired 
walking direction. Related phenomena can not only be 
observed in pedestrian dynamics but also in other phys- 
ical systems. For instance in granular media, the com- 
ponents assort according to their size when vibrations 
are applied to the system @|. The effect occurs particu- 
larly in driven systems Q, e.g., in complex plasmas 
molecular ions [5] and colloidal suspensions [6-8]. For 
the latter also an order parameter for the detection of 
lanes was introduced which we are going to adopt for 
pedestrian crowds (cf. Section III B|) . More on the phase 
diagram and the lane transition of colloids can be found 
inRef. [I[l0|. 

Nevertheless, the phenomenon is best known in pedes- 
trian counterflow from everyday experience although ex- 
perimental or empirical studies are quite rare [ 1 ll — 1 1 5j j . 
The effect can be understood as a self-organization pro- 
cess, where pedestrians try to minimize the contact with 
other pedestrians, especially if they have a different de- 
sired walking direction. The mechanism which leads to 
segregation into distinct lanes is not fully understood. 
A possible explanation could be a fundamental left-right 
asymmetry i.e., people prefer to walk either on the right 
or on the left. But note that this is not necessary for 
the effect. The symmetry can also break spontaneously 
if certain assumptions are made, for instance that people 
tend to follow each other. 

Computer simulations of bidirectional pedestrian 
movement (counterflow) are made frequently 11 5143 ll . 
Lane formation can be reproduced qualitatively [29|, |32j 



and is often used as a validation of a model, but there are 
no quantitative descriptions of the formation process or 
the properties of lanes. Even a description which goes be- 
yond the statement that lanes are present in the system 
or even comparisons with empirical data are extremely 
rare [l5, 33] . Another focus related to bidirectional flow is 
on the jamming transition [UHll], i.e. the transition to a 
state where no movement is possible anymore (gridlock). 
However, there is no empirical evidence for such a transi- 
tion and recent experiments showed that the lower den- 
sity limit for the transition is larger than 3.5 Persons/m 2 
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II. DEFINITIONS 



Floor Field Cellular Automaton Model 



The Floor Field Cellular Automaton (FFCA) model 
[2(il |H, 36] is defined on a two-dimensional square lattice 
where each cell can be occupied by at most one particle 
(pedestrian). In every timestep, each particle is allowed 
to stay at its current position or to move to one of the 
neighboring cells given that the destination cell is not oc- 
cupied by another particle. This is done synchronously 
for all particles (parallel update) . The number of nearest 
neighbors on the lattice can either be four (von Neumann 
neighborhood, cf. Fig.fT]) or eight (Moore neighborhood). 
One can estimate the area assigned to each cell by taking 
the reciprocal of the maximal possible density /) max in a 
pedestrian crowd. Although the empirical value of p ma x 
is not known very well [l[ , a common value often used in 
cellular automata is p max = 6.25 Persons/m 2 which leads 
to a size of 0.4x0.4 m 2 per cell [40l - |42l ]. The length of one 
timestep can be estimated by comparing it with the typi- 
cal reaction time of a pedestrian. This leads to a value of 
roughly 0.3 seconds for one timestep. The movement of 
particles is defined by transition probabilities pij which 
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FIG. 1. Illustration of the transition probabilities in the von 
Neumann neighborhood for a particle located at In 
the Moore neighborhood the four corner cells would have a 
transition probability different from zero. 



(for the model variant [lij studied here) are given by 
Qij 



Pij 



= exp (ksSij + kj^Dij - k A Atj) 



(1) 
(2) 



where fcg and &d are coupling constants. The factor 
£ij G {0, 1} ensures that movement only takes place at 
allowed sites, i.e., £y = if the cell is not acces- 
sible and Cij — 1 otherwise. Inaccessible cells are wall- 
cells, cells which are occupied by other pedestrians and 
cells which are not in the geometry or the corresponding 
neighborhood of the origin-cell. The matrices S and D 
represent the so-called "floor fields" , where S is the static 
floor field (SFF), D is the dynamic floor field (DFF) and 
A is the anticipation floor field (AFF). The latter is the 
main difference to the standard FFC A [20L |35| which cor- 
responds to the choice k^ = 0. The expression 



fen A: 



k\Aij) 



(3) 



in the exponential function of eqn. ([2]) can be interpreted 
as a time-dependent potential, i.e., movement is prefered 
into the direction of larger Vy . 

Bidirectional flow is modelled by two different species 
of particles which have opposite preferred walking direc- 
tions. One type of particle, which is called "type A" , is 
directed towards the right while "type B v particles are 
directed to the left. In a colloidal suspension these types 
would correspond to oppositely charged particles. 

The underlying geometry is a rectangle of W x L cells 
where W and L correspond to the width and the length 
of the system. A subset of L neighboring cells along 
the walking direction is called a row and W neighboring 
cells perpendicular to the walking direction are called a 
column. The particle positions can be described by a 
L x M^-matrix which can be defined as 

-1 if there is a particle of type A 
i (t) = ^ — 1 if there is a particle of type B 
if there is no particle 



(4) 



on the cell The absolute value |r,-j| corresponds to 

the usual binary occupation number. 



1. Static Floor Field 

The Static Floor Field (SFF) does not change in time 
and is not influenced by the presence of particles. It 
contains the information about the preferred direction 
of motion. Typically it describes the shortest distance 
to a destination using some metric. In a countcrfiow 
situation one can give a very simple expression for the 
SFF in equation @ which is 



-1 if the cell is to the right 
-1 if the cell is to the left 
otherwise 



(5) 



for particles of type A and = —S^ for particles of 
type B. This means that different types of particles are 
affected by different SFFs. 



2. Dynamic Floor Field 

The Dynamic Floor Field (DFF) was inspired by the 
motion of ants, who leave a pheromone trace [HI]. Other 
ants are able to smell this trail and follow it. This con- 
cept is adopted in this model: Each particle which moves 
from one cell to another leaves a virtual trace, i.e., the 
DFF value Dij of the origin cell increases by 1. This 
trace acts attracting on other particles due to larger tran- 
sition probabilities according to equation @. The effect 
is that particles have a tendency to follow each other. 
An important detail is to avoid particles being attracted 
by their own virtual trace. Therefore, in the calculation 
of the transition probability in equation ([2]) the DFF is 
temporarily reduced by 1 on the last visited cell. In the 
same manner as the SFF one should use two different 
DFFs Dfj and Dfj which interact only with particles of 
type A and B, respectively. 

The DFF has its own dynamics, namely diffusion and 
decay. After each timestep, it is updated according to 



D iS (t + 1) = (1 - S) D l3 {t) + -AD ZJ (t) 



(6) 



where 



+ D i+1J (t) + A-ij(t) - 4Aj(t) (7) 

is a discretization of the Laplace operator. Equation © 
can also be obtained by discretizing 

d 

g- t D xy (t) = pV 2 D xy (t) - 8D xy {t) (8) 

which is the diffusion equation with diffusion constant 
P = a(l — 8)/ 4 and an extra term for the decay. Here 
V 2 = d 2 /dx 2 + d 2 /dy 2 is the usual Laplace operator in 
two dimensions. 
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3. Anticipation Floor Field 

When pedestrians are in a counterflow they usually try 
to avoid collisions by estimating the prospective route of 
pedestrians with opposite walking direction. This behav- 
ior can be imitated by introducing the anticipation floor 
field [lH • For particles of type A it is defined as 



Ai 



Si T .A dAU ' j,) 



(9) 



where 5%j is the Kronecker delta and A € (0, 1) a param- 
eter which controls the range of anticipation. d,A(j,j') 
is the minimal number of cells which have to be passed 
if a particle of type A goes from a cell (i,j r ) to the cell 
only taking steps towards its desired walking di- 
rection. The field Af^ for particles of type B is defined 
analogously. This definition of Aij causes a large value 
of A* (X 6 {^4, B}) if there are many particles of type 
X which are going to tread the cell in the nearby 
future. Particles of the other type Y should avoid this 
cell, i.e., the AFF acts repulsive (hence the minus sign 
in eqn. @). Note also that only type- A particles are 
affected by the AFF generated by type-B particles and 
vice versa. 



4- Update Rules 

The update rules which describe the transfer from one 
timestep to the next are given as follows: 

1. Computing the AFF: The AFF is computed in the 

beginning of every timestep according to eqn. © 
for both types of particles. 

2. Choosing destination cells: Each particle chooses 

a destination cell according to the transition prob- 
abilities given in eqn. (TTJ and ©. 

3. Solving conflicts: Conflicts are situations, where 

n > 2 particles have chosen the same destination 
cell. In this case one of the particles is chosen at 
random to move. For simplicity it is assumed that 
each particle is chosen with equal probability 1/n. 

4. Movement: Each particle moves to its destination 

cell. If the destination cell is different from the 
origin cell, the DFF value of the latter is increased 
by 1. 

5. Diffusion and decay: The values of the DFF are 

updated according to eqn. (jB]). If the boundary 
conditions are not periodic, it is assumed that the 
DFF is zero beyond the boundary, i.e., Dy = if 
ii {!,■■■ ,W} or j {I,--- ,L}. 

Note that rule 3 for dealing with conflicts can be gen- 
eralized by introducing a friction parameter [20L [35ll36j 



or by coupling with evolutionary games |37H39j . How- 
ever, here we will only study the simplest version as de- 
scribed above since no indication was found that the way 
of dealing with conflicts has a qualitative effect on lane 
formation [44[. 



B. Order Parameter for Lane Formation 

An order parameter which indicates the presence and 
absence of lanes has to detect inhomogeneities parallel 
to the preferred walking directions of the particles. A 
possible approach was given by Yamori [45] who intro- 
duced a band index which is basically the ratio of mov- 
ing pedestrians in lanes to their total number. Other 
attempts for lane detection were made by means of the 
velocity profile [20] and cluster analysis [Hj]. Here we 
will use the same order parameter which has been al- 
ready used in Q to detect lanes in a colloidal suspen- 
sion. The particles there are charged and driven by an 
external field which determines so to speak the walking 
direction. They count for each particle the number JVl 
and No of like and oppositely charged particles whose 
lateral distances, i.e., the projection of distance onto a 
line perpendicular to the field, is below some threshold 
z c (they chose | of the particles diameter). Then the 
quantity cf> n = (A L — No) 2 /(N^ + No) 2 is assigned to 
each particle. This value is close to zero if there is a 
homogeneous mixture of positive and negative particles 
and it is equal to one if there is only one type of particle. 
The global order parameter is given by the average (4> n ) 
over particles. 

This concept can be adopted very well in the FFCA. 
Taking the spatial discretness of the FFCA into account 
we chose z c to be smaller than the diameter of one parti- 
cle. Hence, the calculation of 4> n only considers particles 
in the same row due to the discreteness of space. More 
precisely, to each particle n the value 



m 



N 



N- 



N 



(10) 



is assigned, where i n is the vertical position of particle 
n and Nf (Nf ) denotes the number of particles of type 
A (B) at row i. The global order parameter $ is then 
defined as 



1 N 
N ^ 

n=l 



(11) 



where N is the total number of particles. 

Note that the value of <& is in general larger than zero, 
even if all particles are distributed at random. The mean 
value of $ in that case is denoted by $o- For small den- 
sities $o can be large although there is no lane structure 
in the system. This is taken into account by defining a 
reduced order parameter $ as 

<j> _ (j> n 

$ = — — . (12) 
l-$o 



Note that this order parameter can in principle become 
negative! 

III. DETERMINATION OF STATES 

If not mentioned otherwise, the number of particles 
of type A and B is equal, i.e., pa = Pb = p/2- The 
parameters are given as follows: 

fe S = 2.5, fc D =0, k A = 0, 
a = 0.3, 6 = 0.1, A = 0.8, (13) 
W = 10, L = 100. 

Considering the cell size of 40 cm x 40 cm, this corre- 
sponds to a corridor of 40 m X 4 m. The boundary con- 
ditions are periodic in walking direction and open per- 
pendicular to it, i.e., there is a virtual wall at row and 
PF+1. We use the von Neumann neighborhood, i.e., each 
cell has four nearest neighbors. 
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FIG. 2. (Color online) Definition of states: a) Free flow, b) 
Disorder, c) Lanes, d) Gridlock. 



1. The flow J averaged over the last 50 timesteps is 
lower than 2 ^ L = 5xl0~ 4 . This criterion indicates 
the formation of a gridlock. 



A. Definition of Different States 

The patterns observed in the system can be classified 
into at least four different states, which are defined as 
follows (see also Fig. [5]): 

Free flow: The particle density is very low such that the 
average distance is too large for a considerable in- 
teraction between particles. This leads to a small 
reduced order parameter <I> and a large mean veloc- 
ity. 

Disorder: This state is characterized by a large homo- 
geneity of particles which leads to a small order 
parameter. 

Lanes: Almost every row of cells contains only one type 
of particle apart from a small number of exceptions 
caused by fluctuations. $ is large in this state. 

Gridlock: This denotes a complete jam in the system. 
Particles of type A and B clog each other such that 
no movement is possible in the desired direction. 
The reduced order parameter is usually negative 
and the velocity v is zero. 



B. Simulation Procedure 

The simulations were done as follows: Initially, the 
N = pWL particles are distributed randomly on the 
lattice. Then the simulation starts for different densi- 
ties p e {0.05, 0.06, 0.07, . . . , 1.00} and for different cou- 
pling constants k^ e {0.0, 0.1, 0.2, . . . , 7.5}. For each pair 
(p, ki)) of parameters the simulation was repeated at least 
100 times. The simulation runs were stopped if one of the 
following conditions is fulfilled: 



2. The maximum and minimum value $ max and <I> m i n 
of the order parameter $ in the last 1000 timesteps 
fulfills the following condition 



max ^min 



< o.i, 



(14) 



which indicates the formation of lanes. 



3. The number of timesteps exceeds the density de- 
pendent value T m ax = 20000^/p. This typically 
happens for disordered or free flow states. 

If a simulation is stopped due to the second or third con- 
dition, the velocity v, flow J and order parameter <E> are 
averaged over the last 1000 timesteps. Since the first 
condition means that the average number of particles 
which moved towards their desired direction is less than 
i , it indicates the occurrence of a gridlock. For a given 
set of parameters, the fraction of simulation runs which 
are aborted due to the first condition defines the jam 
probability Pj am - Condition 2 indicates the occurrence 
of lanes. However, it does not contain any restrictions 
for the absolute value of the order parameter because 
more than 99% of all values for $ are anyway larger than 
0.85 given that $ converges. 

Finally, the last condition ensures that the computa- 
tion time for a simulation run is finite. The value of Tmax 
is based on the experience that lanes and gridlocks are 
formed much more quickly, if at all. The factor ^fp takes 
into account that the system takes more time to evolve 
into a stationary state if the number of particles is in- 
creased. Using y/~p instead of p leads to a fast decrease in 
simulation time at low densities where the system does 
evolve neither into a jammed state nor a lane state. 
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C. Results 

1. Jam Probability 

As a first result one observes that in general the jam 
probability -Pjam increases both with an increasing den- 
sity and a decreasing coupling constant ku and re- 
spectively (Fig. [3]). 

This result is not very surprising since an increas- 
ing density means a larger number of opposing particles 
which can block each other. Additionally, there is less 
space to avoid such collisions. An increasing k^ means 
that it becomes more likely for a particle to follow an- 
other particle of the same kind and thus collisions of 
type- A and type-P particles are avoided. That the AFF 
prevents collisions is more obvious, because it prevents 
by definition particles from coming too close to opposing 
particles in the same row. 

The simulations clearly indicate that the AFF is more 
capable to prevent jams compared to the DFF. If the 
AFF is turned off (fc^ = 0, Fig. [3] left) all simulations 
with a density larger than 0.5 evolve into a gridlock. 
Even for p = 0.2 one has a non- vanishing jam proba- 
bility. Anticipation leads to a vanishing jam probability 
for a sufficiently large coupling constant kx > 3 at almost 
every density. 



2. Order Parameter 

Now we consider only those simulation runs which are 
not aborted due to the detection of a gridlock. In gen- 
eral, the averaged reduced order parameter $ becomes 
larger for increasing densities p and coupling constant 
(see Fig. HI). 

As expected, for &d = = the reduced order pa- 
rameter is close to zero for all densities. The fact that it 
is slightly above zero can be explained by the increased 
number of sidesteps a particle has to perform if it has con- 
tact to a particle of different kind. These collisions hap- 
pen more frequently if the number of type- A and type- 
B particles in a specific row is similar. The resulting 
sidesteps lead to a decreasing stability of homogeneous 
states. 

By increasing the coupling constants the order param- 
eter $ shows a specific behavior which depends on the 
density. The dependence is qualitatively the same for 
lanes that are formed by the DFF and AFF, but for the 
latter there are more data for larger densities available 
due to a smaller jam probability. One can see that a 
certain minimal density is needed to reach an order pa- 
rameter which is close to 1. The transition from a disor- 
dered state to a state with almost perfect lanes becomes 
sharper for increasing density. For densities p > 0.2 the 
transition can only be observed at a few data points be- 
cause most of the simulations evolve into a gridlock for 
small coupling constants. In particular there is no stable 



disordered state for large p and the system evolves either 
in a jammed state or in a state with almost perfect lanes. 



3. Dependence on the System Size 

Since the behavior is quite sensitive to the size of the 
system, we want to discuss briefly the influence of the di- 
mensions. The results for the order parameter $, which 
can be found in Fig. show that increases with in- 
creasing width W of the corridor if the length L remains 
constant. This increase seems to converge at a width of 
about W = 30, since the value of the order parameter 
does not change if W is increased further. 

A similar scaling of the system length L could not be 
performed due to the increasing jam probability in that 
case (cf. Fig. As a rough approximation one can 
say that Pj am does only depend on the ratio L/W and is 
monotonously increasing with this ratio. But this is only 
justified for smaller systems and L/W > 1. For large 
values of W one can observe that Pjam is non-monotonic 
whereas Pjam increases with increasing W. 



IV. PROPERTIES OF THE LANE STATE 
A. Distribution of Densities 

In this Section we discuss how the particles distribute 
in the system after the formation of lanes. The data are 
taken after a simulation run was stopped due the lane 
condition (|14[) as described in Section lillBI 

Figure [7] shows the distribution of densities pi in the 
different rows i £ {1, . . . , W}. If the particles were 
equally distributed over the system, one would expect 
a sharp peak in the distribution at the position of the 
global density p. This is for instance the case in a system 
without DFF, even if anticipation is included. But since 
the DFF implies an attractive interaction between parti- 
cles, the number of rows which are not occupied as well as 
the number of rows with larger density grows. This can 
be observed in the distribution as a maximum at pi = 
and a larger support for increasing coupling constant k^>. 
The latter means that the mobility of particles and thus 
the average velocity is reduced. This will be analyzed in 
more detail in Section IIV Bl 

Another feature of the distribution is the existence of 
characteristic lane densities which do neither depend on 
the coupling constant nor on the global density. This 
surprising result was tested for different coupling con- 
stants fcn € [1,7] and for different (global) densities 
p G [0.1,0.3]. These lane densities become visible as 
peaks in the distribution at p\ — 0.28, p\ = 0.46 and 
P3 = 0.56. The different peaks correspond to states 
with different number of lanes, i.e., p* corresponds to 
a configuration with j + 1 lanes. Narrow lanes have a 
higher density than wide lanes. Note that the tendency 
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FIG. 3. (Color online) Jam probability Pj am as function of the density p and the coupling constants A; a an d ko- Left: ko > 0, 
k A = 0. Right: ko = 0, k A > 0. 




FIG. 4. (Color online) Reduced order parameter $ as function of the density p and the coupling constants. Left: ko > 0, 
k A = 0. Right: ko = 0, k A > 0. 

of forming more lanes (-/Vi ancs > 2) increases with increas- type and the maximal possible density in one row of the 

ing coupling constant fcrj. Therefore, for small coupling lane state is pi = 0.5. Therefore, if &d ist large, there is a 

constants (e.g., fco = 1) only 2-lane configurations oc- peak at p* = 0.5 in the distribution. A similar effect can 

cur resulting in a single peak whereas for large coupling be found for large global densities. Then there can be an 

constants (e.g., &d = 5) also 3- and 4-lane configurations additional maximum beyond p% which can even have a 

occur resulting in three peaks. larger value. 

Furthermore, under certain circumstances there are 
additional peaks in the distribution. At small densities 

and large coupling constants fcn almost all particles of B- Velocities 
equal type are located in the same row. This causes ad- 
ditional density dependent peaks. For example at p = 0.1 Since lane formation is a way of separating different 
and a system of size 100 x 10 there are 50 particles of each types of particles which leads to less opposing traffic one 
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FIG. 5. (Color online) Dependence of the order parameter on 
the system size at density p — 0.15. 
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FIG. 6. (Color online) Dependence of the jam probability on 
the system size (p = 0.2, ko = 2). 



would expect that lanes increase the mean velocity v in 
the system, i.e., that v(kA) is monotonously increasing. 
However, this is only true for lanes which are formed due 
to the AFF, i.e., v(kA) is monotonously increasing 44j. 

As one can see in Fig. [SJ the dependence of v on fco 
is not always monotonous. The DFF acts attracting on 
particles of the same type and thus the local densities in 
a row can be large, as seen in the previous section (cf. 
Fig. [7]). This means in general a decrease of velocity, 
because particles are hindered by particles of the same 
kind. Thus there is a decrease in v(kn). Note that a 
state with maximal velocity has not necessarily a maxi- 
mal order parameter. 
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FIG. 7. (Color online) Distribution of densities in the different 
rows. The global density is p — 0.15, the bin size is 0.01. 
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FIG. 8. (Color online) Velocities in dependence of the cou- 
pling constant ko- The data are gained by the simulation 
procedure described in Section Till Bl 



Note that not for all densities a maximum is present. 
For instance at p = 0.1 the global density does not suffice 
to create a local density which is larger than 0.5. Hence, 
each particle has on average at least one empty cell in 
front and the velocity does not decrease. If the density is 
too large one has also no maximum because there are no 
data available for small ko due to a large jam probability, 
i.e., the curve v(kr>) starts at some fcn > and can be 
monotonously decreasing in this case. 



C. Stability and Lifetime 

In agreement with the observations in empirical stud- 
ies [34j] the general structure of lanes in the simulations 
fluctuates in time. The number of lanes changes due the 
vanishing of lanes or merging of neighbouring lanes to 
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a broader one. In this section we will analyze this phe- 
nomenon. Previously, simulations were initialized in a 
disordered state and evolved into a lane state. Now we 
take the inverse approach, i.e., the simulations are started 
in a state where lanes already have formed and we wait 
until this structure is destroyed. The lifetime of lanes will 
serve as a measurement for their stability. Since the sta- 
bility depends on the lane width we use different initial 
conditions with 2, 3 and 4 lanes. The simulation starts 
with one of these initial conditions and is stopped if cither 
the number of timesteps exceeds T max = 1.5 X 10 5 or if 
the order parameter $ falls below . The first condition 
is necessary to perform the simulation runs in an accept- 
able time. This might distort the result for states which 
have lifetimes of order T max . However, T max = 1.5 x 10 D 
corresponds to more than 10 hours in real time since one 
timestep corresponds to about 0.25 seconds. 
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FIG. 9. (Color online) Lifetime of lanes for different initial 
conditions for a density of p = 0.25. 

The results are shown in Fig. GO The two-lane life- 
time increases roughly exponentially with increasing k^. 
It is not surprising that the stability is decreased if the 
number of lanes is increased because then the fraction of 
particles which are at the interface between two lanes be- 
comes larger. The three-lane lifetime only increases up to 
a value of k^> ~ 1.25 where one can see a kink in the curve. 
For larger values of ku, the lifetime stays almost constant 
at a value of about 1000 timesteps. Finally, the four-lane 
configuration has a lifetime which smoothly converges to 
a value of about 250. 

Simulation with anticipation but without the DFF 
show a different behavior. The influence of the number 
of lanes on the lifetime is smaller. The two-lane config- 
uration has still the largest lifetimes, but the difference 
in the stability of the three- and four-lane configuration 
becomes very small for increasing coupling constants. In 
all three cases the lifetime grows faster than exponential 
with Ua- 

The fate of a system which has lost its stability depends 
on the density and the value of the coupling constant 



(cf. Fig. [3] and [4] from Section IIII[) . For large densities 
and a small coupling constant, the system evolves into a 
jam state. For small densities and a small coupling con- 
stant, the system evolves into a disordered state. If the 
coupling constant is large, i.e., at least large enough to 
prevent jams, the system can form lanes again. If the ini- 
tial configuration consists of more than two lanes one can 
sometimes observe that the lanes are merging. In that 
case the final configuration consists of two lanes. Usu- 
ally this happens only in systems with large coupling to 
the DFF, because then the attraction between two lanes 
with the same walking direction is strong. In most cases 
the order parameter $ falls below , such that the sim- 
ulation is aborted and the correct lifetime is taken. Any- 
way, sometimes it happens that lanes are merging with- 
out an appreciable decrease of $ and thus the condition 
for stability is useless. The frequency of these incidents 
increases with increasing k]j and make up to 50%. For- 
tunately, one can easily identify these cases in the data, 
because the lifetime of a 2-lane configuration is much 
larger and thus their is a huge gap (at least one mag- 
nitude) in the distribution of lifetimes. Hence, we only 
take into account the results below this gap. 



V. OPEN BOUNDARY CONDITIONS 

A. Two Groups Passing Each Other 

For this scenario we start the simulation with two 
groups of particles at the very left and the very right 
end of the system (cf. Fig. HUk). Each group contains 
100 individuals of type A and B, respectively. They 
walk in their desired walking direction until they reach 
the boundary of the system. There they are absorbed 
and removed from the system. The simulation ends 
when the flow becomes zero, i.e., if either all particles 
are removed or a gridlock is formed. The simulation 
was repeated 5000 times for each pair of parameters 
(k D ,k A ) G {0,1,2,3,4} 2 . 

The first question is whether the groups are able to 
pass each other or not. If both dynamic floor field and 
anticipation is turned off, the system always evolves in a 
gridlock state. If ko becomes larger the jam probability 
is much reduced, e.g., for kjj = 2 (but kA = 0) about 
50% of simulation runs form a gridlock. For ko > and 
kA > 2 no gridlocks were observed at all. As before we 
only take into account those simulation runs which did 
not end in a gridlock. 

In order to quantify the tendency of collisions, we in- 
troduce a collision index n c as follows: A "collision" is 
defined as a pair of two neighboring cells, where the left 
one is occupied by a type- A particle and the right one by 
a type-B particle. Let iVcol be the number of those pairs 
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in the system. Then the collision index is defined as 

n c = 2— 
where N is the total number of particles. 



(15) 



b) 





d) 



FIG. 10. (Color online) Two groups walking past each other, 
a) Initial conditions. b)-d) Snapshot of the simulation at 
timestep 75. b) ko — 2, k A — 0, c) ko — 0,k A = 2, d) 
kD — kA — 2. 

One can observe that the dynamics shows qualitative 
differences depending on the choice of the coupling con- 
stants. Without anticipation it is strongly influenced 
by collisions between opposing particles (cf. Fig. [TOb 
and [TT|) . Even after most of the particles have already left 
the system there are almost always small groups of par- 
ticles which form a jam that dissipates after longer time. 
Therefore, the collision index n c fades very slowly. Lanes 
can only be formed after both groups get in direct contact 
with each other, but they are less distinctive compared to 
those in the periodic system. This leads to a quite small 
order parameter. With anticipation collisions are much 
less likely, but there are also difference in the dynamics 
depending on whether the DFF is turned on. If not, many 
thin lanes are formed during the contact fFig. [TUbh Af- 
ter the groups have passed each other these lanes are de- 
stroyed immediately. This can also be observed in terms 
of the order parameter in Fig. QTJ Probably the most 
realistic results were attained by a combination of AFF 
and DFF. The tendency of collisions is reduced further. 
Lanes can be seen clearly even if the order parameter is a 
bit reduced compared to the system with AFF only. This 
can be explained by lanes which are often not parallel to 
the walking direction fFig. [TUtO. Their number is much 
decreased and one usually observes 2 or 3 lanes. The lane 
structure is kept after the groups have passed each other. 
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FIG. 11. (Color online) Order parameter $ and collision index 
n c from eqn. (|15[) versus time for two groups walking past each 
other. 



only type-B particles whereas the right boundary absorbs 
only type- A particles. 
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FIG. 12. (Color online) Jam probability at timestep 20000 
for open boundary conditions. 



B. Random Insertion of Particles 

Now the particles are inserted and removed dynami- 
cally with a certain insertion probability f3- m — /3 on each 
cell next to the boundaries. At the left (right) bound- 
ary only type- A (type-B) particles are inserted. Like in 
the previous section, the removal of particles will always 
happen, i.e., with probability /3 out = 1 at the left and 
right boundary, respectively. The left boundary absorbs 



Fig. rT2] shows the probability that the system evolves 
into a gridlock before timestep 20000. It turns out that 
the avoidance of jams is much harder as in the periodic 
system. The DFF does not help very much to prevent 
jams but astonishingly it rather supports the formation 
of a jam. 

Anticipation helps much to avoid gridlocks and even 
for an insertion probability of /3 = 0.12 gridlocks can be 
avoided. We find a linear dependence between f3 and p. 
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The constant of proportionality depends slightly on the 
exact value of the coupling constants and ranges from 
2.5 to 2.8. This means that the maximal density in the 
system without forming a gridlock is about p — 0.34. 

Before gridlocks are formed, the system usually forms 
lanes for sufficiently large coupling constants. The order 
parameter is lower than in the periodic case, but compa- 
rable to the findings in the previous section (Fig. ITT]) . 

It is not surprising that systems with open boundary 
conditions show a slightly different gridlock behaviour 
from that with periodic boundaries. In the periodic case 
the flow in a system near a gridlock becomes very small. 
In contrast, in the open system still additional parti- 
cles are fed into the system. Due to particles moving 
up " from behind" the chances for the gridlock to resolve 
are strongly reduced. 

VI. CONCLUSIONS 

We have studied the phenomenon of lane formation 
in a cellular automaton model for pedestrian dynamics. 
In contrast to previous studies, we have obtained for the 
first time quantitative results. This was possible by intro- 
ducing an order parameter which has been adapted from 
the analysis of colloidal suspensions. It has been shown 
that this order parameter is suitable for the detection of 
lanes, at least for the simple scenario of a straight corri- 



dor considered here. 

It is found that at least four different states can be 
distinguished, namely free flow, disordered flow, lane for- 
mation and gridlock. Using the order parameter we have 
mapped out the parameter range of the model for which 
lane formation is possible. Also other quantities like den- 
sity and velocity show a very characteristic behavior in 
these scenarios. 

Most models for pedestrian dynamics show a strong 
tendency towards gridlock which appear for large densi- 
ties. This is rather unrealistic and was not observed in 
experimental studies. A model extension which includes 
an anticipation mechanism was found to show a more re- 
alistic behavior by suppressing the formation of gridlocks, 
at least in the case of periodic boundary conditions. 

Although we have focussed on the floor field model and 
its variants we want to emphasized that the techniques 
developed here can also be applied to most other models. 
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